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^ . Abstract 

I Entanglement types of pure states of three spin-i particles are 

' classified by means of their stabilisers in the group of local unitary 

^ transformations. It is shown that the stabiliser is generically discrete, 

. . i and that a larger stabiliser indicates a stationary value for some lo- 

I cal invariant. We describe all the exceptional states with enlarged 

^ ' stabilisers. 

1 Introduction 

It is only relatively recently that the importance of entanglement has been 
fully realised. Not only, as Schrodinger emphasised does it constitute 
one of the chief differences between classical and quantum mechanics, and 
the main obstacle to an intuitive understanding of quantum mechanics; the 
recent discovery is that it is also a resource, yielding much greater capabilities 
than classical physics in information processing and communication (see for 
example 0). 

It is therefore important to analyse and measure this resource. A full anal- 
ysis has so far been achieved only for pure state systems with two component 
parts P, H]; for multipartite systems there are several different possible mea- 



sures of entanglement [pi 0, 0, U, 00, 11 1, the relation between them being 
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incompletely understood. A full quantitative analysis of entanglement even 
for pure states of three-part systems appears to be difficult (but see 0). 
Our aim in this paper is to give a qualitative analysis of the entanglement of 
such states, using group-theoretic methods to classify the possible kinds of 
entanglement. 

The nature of the entanglement between the parts of a composite sys- 
tem should not depend on the labelling of the basis states of each of the 
part-systems; it is therefore invariant under unitary transformations of the 
individual state spaces. Such transformations are referred to as local uni- 
tary transformations, though there is no implication that the part-systems 
should be spatially separated. If the part-systems have individual state 
spaces Hi, . . . , 7i„, so that the space of pure states of the composite system 
is Hi®. . .^TCn, then a local unitary transformation is of the form Ui®...®Un 
where f/, is a unitary operator on Tii. The set of all such transformations is 
a group G, whose orbits in Hi ® ... ® 'Hn are equivalence classes of states 
with the same entanglement properties. Each orbit therefore corresponds 
to a complete specification of entanglement. The orbits can be classified by 
their dimensions, which are determined by the stabiliser subgroups of points 
on the orbit; the relation is 

dimO + dim>S = dimG (1) 

where O is an orbit and S is the stabiliser of any point on O, i.e. the set of 
elements of G which leave a point unchanged (different points on the same 
orbit have conjugate stabilisers, which have the same dimension). 

This paper is concerned with pure states of three spin-| particles (n = 3; 
Til = ^2 = Hs = C^.) We will show that for most states (all but a set of lower 
dimension) the stabiliser is discrete, so the dimension of the orbit is the same 
as that of the group G. Classifying types of entanglement by the dimension 
of the orbit is therefore equivalent to identifying certain exceptional types 
of entanglement, which can be expected to be particularly interesting and 
important. One way in which this manifests itself is that any such exceptional 
entanglement is necessarily associated with an extreme value of one of the 
local invariants which form coordinates in the space of entanglement types, 
and from which any measure of entanglement must be constructed. 

The organisation of the paper is as follows. In Section 2 we review the 
case of two spin-i particles. The results here are well-known, but we include 
them for the sake of completeness and orientation. In Section 3 we prove 
the general theorem about three spin-i particles mentioned in the preceding 
paragraph. Section 4 consists of the theorem concerning the association 
between enlarged stabilisers and stationary values of invariants. Section 5 
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contains the classification of exceptional entanglement types in the system of 
three spin-i particles, in which we examine all the states which are identified 
as non-generic in the theorem of Section 3. Section 6 is a summary listing 
these exceptional states. They are illustrated by means of plots of their 
two-particle entanglement entropies in an appendix. 
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2 The Stabiliser for the 2-particle case. 

A pure state of two spin-^ particles can be written as 

\^) = J2U,\^P,)\^|J,) (2) 

i 

where IV'2)} is a basis of one-particle states. Having fixed this basis, 

we can identify the state with the matrix of coefficients T = (tij). The 
group of local transformations is 

G2 = U{1) X SU{2) X SU{2), (3) 

since the phases in the individual unitary transformations can be collected 
together. The effect of a local transformation (e*^. A, F) on T is to change 
it to e'^^XTY"^, so the condition for (e*^. A, F) to belong to the stabiliser of 
1^') is 

T = e'^XTY^. (4) 

For a 2-particle state, we can always perform a Schmidt decomposition, so 
we need only consider states for which 




i.e. 



^)=p| T)IT) + g| i)U) 
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where p, q are real and positive. Multiplying the stabiliser equation on the 
right by Y, where the overbar denotes complex conjugation, and writing 

we obtain: 

fp 0\f f ( g_ h\ fp 0\ 

\0 qj \-s r) \-h g) \ Q q) ' 

For given p, q, we want to find the set of solutions (</?, h, r, s) with (p real 
and g, h,r,s & C, with \g\'^ + \h\'^ — 1, |rP+ |sp = 1. If p 7^ 0, then g — fe"'^'^ . 
If g 0, then 'g — re~"^. Therefore either r = or (/? = wk. Also 

h = ^se-''^ = -se''^ 
q p 

So unless p — q (since p and q were obtained by a Schmidt decomposition, 
they cannot be negative) we must have ^e"**^ — ^e*"^ 7^ and so s = 0. The 
states now fall naturally into three classes: 

Case 1: The General case. If p 7^ and g 7^ and p ^ q then s = /i = 
and e**^ = e~'^f = ±1 so we can absorb that external sign into X. This is the 
subgroup 

(^ = 0, X = e''"'% Y^X. (6) 

The stabiliser has one parameter, v. 

Case 2: The Unentangled case. Without loss of generality, we can take 
p = 1, g = 0. Putting g = e^^, this is the subgroup 

(e^^^,r) = (e^^e^^e-(^+^)) (7) 

The stabiliser has two parameters, ip and 9. 

Case 3: The Mciximally Entangled Case. This occurs when p = q — 
I/V2. Then 

g = re^*'^ = re*"^ 
h = -se-'^ = -se'^P 
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So if = UTT (or else we'd have to have r = s = which is impossible). 
Thus g = ±r and h = =ps, giving the three-parameter subgroup defined by 
Y = ±X, where X can be anything in SU{2). 

These results illustrate how the occurrence of a state with special phys- 
ical significance is signalled by a change in the stabiliser. In Case 2 above 
the states are factorisable, so there is minimal entanglement: the stabiliser 
increases from one- to two-dimensional. In Case 3, on the other hand, the 
entanglement is maximal as measured by the entropy of entanglement 

S = \np^ + In 

or equivalently by the 2-tangle |Q 

2 2 2/1 2\ 

T = p q = p [1 — p ) 

(see Section 6). We note that the stabiliser for these states is even larger, 
being three-dimensional. 

This association between an enlarged stabiliser and a maximum or min- 
imum of an invariant measure of entanglement is a general phenomenon, as 
will be proved in Section 4. 

3 The 3 spin-^ Particle Generic Stabiliser. 

In this section we will show that the generic pure state of three spin-| par- 
ticles has a discrete stabiliser in the group 

G3 = f/(l) X SU{2) X SU{2) X SU{2) (8) 

of local unitary transformations. This is in contrast to the case of two par- 
ticles, where, as shown in the previous section, every state has a stabiliser 
which is at least one-dimensional. In the course of the proof we will iden- 
tify those exceptional states for which the stabiliser might have dimension 
greater than zero. For ease of later reference, we will label those steps in the 
argument whose failure could produce such nongeneric behaviour. 

Theorem 1. Let |\E') be a pure state of three spin-^ particles, and let L(^) 
be the Lie algebra of the stabiliser of |\E') in the group G3. Except for a set of 
states |\E') whose dimension is less than that of the full space of states, 

L(^) = 0. (9) 
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Proof. Any state of three spin-| particles is of the form: 

\'^)-Yl^i3k\'^iM)\'tpk) 

where i,j,k = 1 or 2 and j'^i) = | T); |'02) = I i)- A local transformation 
is of the form: 

|*)l-^e''^ ^ tijkUeiVmjWnkHe)\lpm)Hn) 

for some 2x2 matrices U,V,W e SU (2) and some phase (p. Suppose U, V, W 
are close to the identity: 

U = 1 + isA, V = 1 + isB, W = l + ieC (10) 

where e is infinitesimal and A,B,C are hermitian and traceless. li 9 — e(p is 
also small we have, to first order in e, 

ijk + (^ietejk + bjmtimk + Ckntijn) I V'i ) I V^j ) I V^fc ) 
= ^ {{^^ii + aiijUjk + bjmtimk + CkJijn) I'ipi) IV'fe) 

Hence if the local transformation (e*'', U, V, W) belongs to the stabiliser of 

{ifiSii + aii)tijk + bjmtimk + Ckntijn = 0, (11) 

using the summation convention on repeated indices. Let Tj be the matrix 
whose (j, fc)th entry is tijk] then these equations can be written in matrix 
form as 

(ifiSu + au)Tt + BTi + TIC^ = 0. (12) 

Separating these at their free indices, and performing the summation gives: 

STi + 71 + ((^ + aii)Ti + 012^2 = 
BT2 + T^C^ + a2iri + (^ + 022)^2 = 0. 

Generically (Gen 1), at least one of Ti and T2 is invertible (say T2); if so, 

BT{r2-^ + TiC^T2-^ + (ip + an)TiT2-^ + au ^ 
TiT^-'B + TiC'^T^-^ + {ip + a22)TiT2-' + 021 (71^2"')' = 

6 



and so 

-TiC^Ta-^ = ETiTs"^ + + aiOTiTs"^ + 

= TiT2-^B + {ip + a22)TiT2-^ + aai (TiTs"')' 

Let X = TiT2~^; then these equations give 

[B, X] = -ai2l - (an - 022) X + a2iX^ 

Now we use 

tr (X" [B, X]) = tr (X'^EX - X^'+^B) = 

to obtain 

tr (012! + (an - 022) X - a2iX^) = 

and 

tr (012^ + (an - a22) X"^ - a2iX^) = 0. 

Let A and /j, be the eigenvalues of X. Then we obtain 

2ai2 + (an - a22) {X + /i) — a2i (A^ + A*^) =0 
ai2 (A + //) + (an - a22) (A^ + ii"^) - a2i (A^ + /x^) = 

Generically (Gen 2), A + 7^ and so solving for ai2 and a2i in 
(an - a22) gives 

Xfi 1 
ai2 = -T— — (an - a22j , a2i = t— — (an - a22j , 

A + H A + H 

but generically (Gen 3) this will not satisfy ai2 = 021 unless 

ai2 = a2i = (an - a22) = 0. 

Hence A — 0, and the equations for B and C become 

{B - (fl) Ti + TiC^ = 0, 
BT2 + T2 (C^ - (pi) = 0. 

The second of these equations gives C as 

= -Tr' {B - ipl) Ti = -T2-^BT2 + ipl. 



Taking the trace of this equation, (p = 0. Putting this into the first equation 
shows that B commutes with TiT2^^. Generically (Gen 4), the only matrices 
which commute with a 2 x 2 matrix X are al + jSX for some scalars a, (3] 
therefore 

B = al + (3T,T2-\ (26) 

Generically (Gen 5), this will not be hermitian unless (3 = 0, and then 
tri? = implies a = 0. Thus B = and therefore C = 0. Thus for generic 
values of tijk the only solution of (|T2|) is 

^l = A = B = C = 0, (27) 

so the stability group is discrete. □ 

Remark: It follows from this theorem that the generic orbit has the same 
dimension as the group G3, namely 10. Since the space of (non-normalised) 
state vectors has (real) dimension 16, the number of independent invariants, 
which is the same as the dimension of the space of orbits, is 6 (including the 
norm) . 



4 Exceptional States: The significance of an 
enlarged stabiliser 

In this section we will prove that a three-qubit state which is exceptional 
in the sense of Theorem 1 has a stationary value of some fundamental in- 
variant. Since any measure of entanglement must be such an invariant, this 
indicates that these mathematically exceptional states are likely to have a 
special physical significance. 

By a "local invariant" we mean a real-valued function of the state vector 
which is invariant under local unitary transformations, and is therefore con- 
stant on each orbit. It is convenient to concentrate on polynomial functions, 
which can be regarded as coordinates on the space of entanglement types; 
more general invariants (e.g. entropy of entanglement) can be constructed 
from these. Since the generic orbit in the state space Ti has dimension 
dimGs, the number of parameters needed to specify such an orbit is dimH— 
dimGs = 6. Such parameters, being constant on orbits, are invariants. 

The space of orbits is not necessarily fiat, and it may not be possible to 
parametrise it globally with a single set of six invariants (see 0): geomet- 
rically, the space of orbits is a manifold which may have several different 
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coordinate patches; algebraically, the algebra of invariants is not a polyno- 
mial algebra but is generated by more than six invariants which are subject 
to some relations. However, we can choose a neighbourhood of a state so that 
the algebra of invariant functions on that neighbourhood has six independent 
generators. 

Theorem 2. Let Ti he the space of 3-qubit pure states, and let G3 be the 
group of local unitary transformations ofTi. Let Ii, be a set of 6 polyno- 
mial invariants which generate the algebra of local invariants in a neighbour- 
hood of a state {ipo)- U the stabiliser of \iPq) in G3 has non-zero dimension, 
there is a linear combination 0//1, Iq which has a stationary value at \ipo). 

Proof. Let xi, xie be real coordinates on H. Suppose the Jacobian matrix 

has maximal rank 6 at lipo). Since the are polynomials, the 6x6 minors of 
J are continuous functions, so if one of them is non-zero at \ipo) it is non-zero 
in a neighbourhood of {ipo)- Hence, by the implicit function theorem, the 
equations 

/,(|^)) = /,(|^o)) (29) 

define a smooth manifold in Ti. of dimension diniTi— 6. These are the equations 
of a level set of the polynomial invariants of G3. Since is compact, its 
invariants separate the orbits [0 and so (p9D is the equation of the orbit 
of {ipo), which therefore has the same dimension as G3. It follows that the 
stabiliser of \ipo) is discrete. 

Hence if the stabiliser of {ipo) is not discrete, then the matrix J has rank 
less than 6 and therefore there exist scalars (Ai, Xq) such that 



i.e., the linear combination 



(31) 

has a stationary value at \iPq). □ 

Note that this theorem does not guarantee that all stationary subspaces 
of any invariant will be associated with enlarged stabilisers. However, it does 
indicate that states with enlarged stabiliser dimensions are likely to have 
special physical significance. 
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5 The classification of non-generic states 



5.1 Setting up the problem. 

We will look for the stabilising subgroup of the group G3 = f/(l) x 5'f/(2)^ 
of local transformations, i.e. the group of (e**^, ?7, V, VT) where U,V,W are 
all elements of SU{2) and e*"^ is an overall phase. We will start with the 
three-index tensor equation for the local transformations: 

t ijk '^^^^'^il'Vjm'^kn^lmn (^2) 

where the t's are the coefficients of the state vector and the uus are the 
matrix elements of f/ G SU (2) etc. Using the (Ti)-^ notation introduced in 
Theorem 1, and partitioning the equation at the index i: 

T\ = e'^PV [uuTi + M12T2] (33) 



T'2 = e'^py [naiTi + U22T2] (34) 

where U22 = uii and U21 = —U12 and + |Mi2p = 1- The stabiliser is the 
set of (e^^, U, V, W) such that T'l = Ti and T'2 = T2. 

In examining the non-generic states, not covered by Theorem ^ whose 
stabilisers have potentially non-zero dimension, we will sometimes find it con- 
venient to abandon the infinitesimal approach of Theorem 1 and determine 
all finite elements of the stabiliser groups. 



5.2 The "bystander" rule. 

We will now examine the apparently trivial case when either Tj (say Ti) is 
the zero matrix. In this instance it is possible to choose bases of the two 
one-particle spaces such that Ti and T2 become diagonal. We need therefore 
only consider the case 



a 

p 



where (3 may or may not be zero. Then the first stabiliser equation (|33|) 
becomes 

= e'^^V (MnO + U12T2) 
therefore uu = U21 = 0, and the other equation becomes 

T2 = e'^^V {e^'^T2) (35) 
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where e'^ = 1^22- This can be seen to be the 2-particle stabihser equation, but 
with an additional external phase factor - which for the sake of transparency 
later we will not absorb into ip. The fact that one of the Tj-matrices is the zero 
matrix means that states of this type are factorisable. The particle(s) whose 
kets can be factored out in this way do not participate in the entanglement 
(if any) of the other particles and so we'll call these 'bystander' particles, 
and states in which not all the particles participate in the entanglement 
'bystander' states. 

If T2 is singular, we have the equation 

« 0\ ^ , fa 0\ T 

which, by the two-particle result reduces toV — e^'^'^^ , W — e*^'^^ ^j^h 

a = e^'^e'^e^^ae^" (36) 

giving us the condition 

i.e., three degrees of freedom. 

If T2 is non-singular, use Section 2 to look up the appropriate 2-particle 
stabiliser. This comes down to whether or not \a\ — If \a\ 7^ the 
equation becomes 



(3 = e^'^e^^(-e^^)/3(-e^'') 



which both reduce to 



ip + 6 + -f + r]^ 2mr, (37) 



which makes three degrees of freedom. 

If \a\ — \P\ , take a = /3 e R. Then we have 



al = e'^e'^^UalU^ (38) 



so 

(f + 9 — 2mr 

and one element of SU(2) giving us four degrees of freedom. 

Thus (in the three spin-| case) factorisable states reproduce the stabil- 
ising group structure of the fewer-particle states that their sub-systems re- 
semble. 
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5.3 Exchanging the particle labels. 

Recall that in Theorem 1 we chose particle 1, with corresponding index i, as 
the 'partitioning index' which splits the original, 3-index state vector 'tensor' 
problem into the more manageable form of a pair of coupled matrix equations. 

'Pi '■ Ujk {Ti)jk (39) 

This choice of particle 1 was entirely arbitrary: we could just as easily have 
chosen either of the indices j or k. Changing the partition index is sometimes 
useful. The effect of changing the particle labels (repartitioning) on the 
stabiliser is simply to permute [/, V, W as each particle's associated SU (2) 
copy just follows its associated index. 

In group theoretical terms, the operations of permuting the particles are 
unitary operations on three-particle states which, though not elements of 
the group of local unitary transformations, do belong to the normaliser of 
this subgroup in the group of all unitary transformations. States related by 
elements of the normaliser will have isomorphic stabilisers in the group of 
local unitary transformations. 

5.4 Change of basis 

We are, of course, always free to change the basis that we use to describe 
states of any of the three particles. (This amounts to applying a local uni- 
tary transformation in the passive interpretation.) If the change of basis is 
described by the 2x2 matrix P for particle 1, Q for particle 2 and R for 
particle 3, then the effect on the matrices [/, V, W is 

U^PUP^\ V-^QVQ^\ W^RWR-\ (40) 

The effect on the matrices Ti,T2 is the same as in (|33D, ([3^ ) with {P,Q,R) 
replacing {U, V, W). In other words, the group element (e*"^, U, V, W) is conju- 
gated by the group element corresponding to (P, Q, R) (namely (e*^, P', Q', R!) 
where e'^ = (det P ■ det Q ■ det R)^ and P' = (det P)'^P, etc.) 

If we regard P,Q,R as active transformations, taking the state (Ti,T2) 
to a different state on the same orbit, then this is the basis of our earlier 
remark that all the points on a given orbit have conjugate stabilisers. 

5.5 Type 1 Non-Generic States: Both T-matrices sin- 
gular 

In Theorem |l] the first step in the argument that is only generically true 
(Gen 1) needs at least one Tj to be invertible for the argument to be valid. 
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If both Tj's are singular, we can choose our coordinates to put one Tj, Ti say, 
into diagonal form by an appropriate local transformation. Then Ti and T2 
will be of the form 

where the singular value p is real and positive (the case where p = has 
already been dealt with in subsection 5.2.) The stabiliser equations, obtained 
from (|33|) and (|3^) by imposing the conditions that T'l = Ti and T'2 = T2, 
are: 

Ti = e'W [uuTi + U12T2] (42) 



T2 = e'W [UnTa - II12T1] W^. (43) 



From (^) and ( ^31) it can be seen that a necessary condition for an 
enlarged stabiliser to occur is that muTi + U12T2 and — II12T1 + M11T2 must 
have the same singular values as Ti and T2 respectively. In particular, they 
must have the same determinant, namely zero. Taking the determinant of 

MllTi + M12T2, 



uuuupbd = 0. (44) 



We will write 



5.5.1 Case 1: a,b,c,d all nonzero (Semigeneric states) 

Suppose a, b, c, d are all non-zero. We will call this form "Semigeneric" , as 



it is the generic form for a singular matrix for T2. Equation (|4^) shows that 
either uu = or U12 = 0. If Uu = 0, write Uu = e*^; then (^) becomes 

P 0\ 5fr -gs 



oy-P»""^"'(-& fej- (^«) 

Hence h = s = 0, (7 = e*", r = e*'^, with 

If + 6 + a + (3 = 01271. (47) 
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Now gives 

\hc hd) ~^ \^e*(-"+^)6c e-'^'^+'^^bd 

so 

a + (3 = a-(3 = -a + (3 = -a-(3 = e-ip (mod 27r). (49) 

From this, together with (0), it follows that each of the angles if, 6, a, (3 is 
equal to or vr and therefore the stabiliser is discrete. 

We will write the stabiliser as 5 = SiU S2, where Si is the subset with 
U12 = and ^2 is the subset with uu = 0. Then Si is a subgroup. The 
product of any two elements of ^2 belongs to Si, so ^2 is a single coset of Si 
(unless it is empty) and therefore contains the same number of elements as 
Si, and is therefore also discrete. 

Case 2: a = or b = 0, bd ^ (Slice states) If either a or c = and 
the other three of a, b, c, d are non-zero, then the state is either 

P\m) + bc\ii]) + bd\iii) (50) 
or p|tTT) + arfUn) + M|iii) 

which are equivalent to each other under exchange of particles 2 and 3. (The 
third similar state, 

p\m) + q\ni)+r\lU) (51) 

can be obtained by a permutation of the particle labels.) For the state (|50| ) 
the equations for uu = give the one-dimensional set of stabiliser elements 

{e'"^, U, V, W) = {eil, e^^"^ £2e"^^"^ £1621) where 61,62 = ±1. (52) 

The equations for uu = require Ti and T2 to have the same singular 
values, the condition for which is 

/ = l&nicP + m. (53) 
If this is satisfied, the stabiliser equations are 
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These give the stabihser elements with un = as 



{e'^,U,V,W) = 



p p 



(55) 



where x =arg(&c) and 6 can take any value between and 2n. 

Thus the slice states have a one- dimensional stabiliser consisting of the 
four circles (|^) unless (|^) is satisfied, when the stabiliser is doubled and 
also contains the four circles (^). We call this set of states a "slice ridge". 



Case 3: a = c = 0, bd (The GHZ states) If a = c = 0, but 6c/ ^ 
the state is the GHZ state. 



p\m) + q\lU) (56) 

with p and q = bd both non-zero. We may assume that they are both real 
and positive. The singular value condition tells us that unless |g| = p the 
only solutions to the stabiliser equations will have Ui2 = 0, giving the two- 
dimensional stabiliser 

{e'^,U,V,W) = (±l,e^^"^e^°"^e^'^"^) 

with the condition that 9 + a + P = 0oYn. 

If |g| = p, the stabiliser is doubled, and also contains the elements 



{e'"^, U, V, W) = 




with the condition that 



9 + a + p = 0oni (57) 

This is the original GHZ state, which can be regarded as a three-particle 
analogue of the maximally entangled ("singlet") two-particle state. We note 
that although the GHZ state has an enlarged stabiliser when its coefficients 
are equal in magnitude, the enlargement does not consist of an increase in 
dimension as in the two-particle case. 
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5.5.2 Case 4: 6 = or = (Bystander states) 

If 6 or or both are zero, the determinant equation (|4^) no longer imphes 
that U must be either diagonal or ant i- diagonal. However, in all of these 
cases the state factorises and one of the particles is a bystander. We will 
just look at the 6 = d = can be obtained by the appropriate 

transpositions, and go back to the "both" case after that. We have the state 
vector: 

_/pO\ rr _ (ac ad\ 

i.e., 

\h){p\ TiT3) + ac| U],) + ad\ ilia)) 

which is a state in which particle 2 is a bystander, and therefore has been 
dealt with in section 5.2 above. 



5.5.3 Case 5: h = d = Q (completely factorised states) 

In this case, 

_/pO\ _ fac 

so the state vector is: 

ip\ h)+ac\ T2T3) (58) 

which is the totally factorised state, and has already been considered as the 
T2 singular bystander case. 



5.6 Non-generic Type 2: tr(rir2"^) = 0. 

Let us now consider what might happen if the assumption (Gen 2) fails. If 
A + /i = 0, equations ([TTD and (|T8|) become 

ai2 = A^a2i, (59) 
2A2(an -022) = 0. (60) 

We can still deduce that A = (since ai2 = 021 and an + 022 = 0) unless 
|A| = 1 or A = = 0. 
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5.6.1 Case 1: |A| = 1. 

Since = 021, equation (p9D gives = aX where a is real. The right hand 
side of ([ISD becomes 

aX{X^-XH) (61) 

by the Cayley-Hamihon theorem. Thus it is still true that B must commute 
with X. We can change basis for particle 2 (multiplying Ti and T2 on the left 
by a unitary matrix P) so that X takes the form 

X = = . (62) 

Since X is not a multiple of the identity, the requirement that B should 
commute with X gives 

B = ul + vX (63) 

for some scalars u, v; but B is traceless, so m = 0. 

Suppose u 0. Since B is hermitian, f = 0; thus B = 0. Now equation 
dH) gives 

= -^i - aATs^^Ti. (64) 

Hence 

^ = -^tr [C^ + «AT2-^ri] (65) 
= -^tr [aXT^T^-'] = 0. (66) 

Now we can change basis for particle 3 (multiplying Ti and T2 on the right 
by a unitary matrix) so that T2 takes the form 

with a 7^ since T2 is invertible. Then 

CT _ -aXT2'^XT2 (68) 

A a-i(Acj + 26)\ 
= -« L M- (69) 
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Since C is hermitian, a non-discrete stabiliser can only occur if 

A^-^. (70) 

Then the state is 

l^) = A| T)(a| T)IT)-&I T)li)-Ii)li)) 

+ U)(«IT)IT) + &IT)U) + U)U)) (71) 
= a|T')IT)T) + &ir)IT)li)-A|nii)|i) (72) 

where 

in = ^(Air) + ii)), 
in = ^(iT)-Ai i)). 

This is one of the slice states considered in Section 5.5. 
If a; = 0, equations (^) and (|67D immediately give 



l^) = A|T)(a| T)|T) + 6| T)U)-U)U)) 

+ li)(«IT)IT) + &IT)li) + li)li)) (73) 

= a|niT)T) + &U')lT)li)-A|r)U)U) (74) 
which is again a slice state. 

5.6.2 Case 2: A = /i = 0. 

The only remaining possibility is that X = TiT2^^ is unitarily equivalent to 

oj- 

In this case (|1^) and ([T8|) give only ai2 = 021 = and (|15D becomes 

BX -XB = 2rX (76) 

where r = an = —022- With X = ' follows that 
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i.e., B = A. Now we return to equations (|T2|) of theorem |1|: 



Writing 



so that 



and 



these become: 



{ipSii + au)T, + BTi + T.C^ = 0. 



a b 
c d 



r 
-r 



ujc ujd 




+ 



and 



r 
-r 



a b 
c d 



+ 



UJC ujd 




a b 
c d 



UJC ud 




s y 

y -s 



X y 

y -X 



X y 

y -X 



-(V3 + an) 



- an) 







a b 
c d 



These give us the following constraints: 

c(s + r) + dy = —c{ip + r) 
d{r — s) + cy = —c{ip + r) 
a(r + s) + by = — a((/? — r) 
b{r — s) + ay = —b{ip — r) 
c{s — r) + dy = —c{ip — r) 
—d{r + s) + cy = —d{ip — r). 

which produce just four independent equations: 

^{f + r) -^if -r) - as - by 

b, N 36, , , 

-{if + r) - —{if -r)+bs-ay 

c c 

-{if + r) + -((^ -r) + cs + dy 
+ r) + ^((^ -r)-ds + cy 








0. 
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For a non-zero solution ((/?, r, s, y) with (p, r, s real, the matrix 

/a -3a -2a -2b\ 

b -36 26 -2a 

c c 2c 2d 

\d d -2d 2c ) 

must have determinant zero. This gives us that 

det(r2)ac + det(7^)6d = 0. 



(82) 



(83) 



Since is non-singular by assumption, this allows us only three possible 
solutions: 



a = (i = 
c = 6 = 
|ac| = |6(i|, all non-zero. 

If a = d = we have by — therefore y — 0. Then 



(84) 
(85) 
(86) 



{r — s) = -{(f — r) 
{s — r) — —{ip — r) 

which gives us that {(p — r) — and also that s — r. This solution has one 
degree of freedom, which we'll call (p. The state is: 



ujc 




To 



6 
c 



and the stabiliser for states of this type is. 

If 6 = c = we have a state vector that looks like this: 



(87) 



(88) 



ujd 




To 



a 
d 



which is just a reflection of the state vector in the previous case in the vertical 
midlines, and so can be mapped into it by a change of basis, as can its siblings 
obtained by permuting the particle labels. The stabiliser for these is thus: 
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so relabelling the spin coordinate just relabels the stabiliser variable, as ex- 
pected. We nickname these states "Beechnut" states, because when the three 
one-particle von Neumann entropies for this subspace are plotted, we think 
it looks like a beech nut. 

This leaves us with the "non-zero" solution. It can be seen that {ip + r)c = 
dy and hence that 2rc = — 2rc which means that r = since we've assumed 
that c 7^ 0. Hence r = s = and by = (fa. So we have 

by = ipa 
_ a 

^=^^ 

and so 

a c 
b^d' 

Therefore 

ad = be 

and the determinant of T2 is zero after all: this case is Type 1 Non-generic, 
and is in fact a bystander case. 



5.7 Non-generic Type 3 

In this next stage of the calculation, we will assume that both Ti and T2 are 
non-singular, and move on to consider the failure of the assumption (Gen 
3). In Theorem 1 we obtained the equations (pi)) 

ai2 = - , (an - 022) 021 = — (an - 022) (89) 

A /i X + fi 

where A,/i are the eigenvalues of the matrix X = TiT2~^. But generically, 
this will not satisfy a^ = 021 unless 

ai2 = 0-21 = o-ii — 022 = 0, 

so that A = 0. We will now examine values of A and fi that allow A to be 
non-zero. 

Since A is hermitian and traceless, an = —022 is real. So ai2 = 021 
requires 

_ A/x _ 1 

A + /i A +71 
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I.e., 

— |A|^/i — A|/ip = A + /i. 

Now we know that \X^\ = 1 from these same equations. Substituting for 
gives 

(|A|2 + 1) (/i|Ap + A) =0 (90) 

Hence 

A (A/i + 1) = (91) 
and so the eigenvalues of TiT2~^ must be of opposite phase, namely: 

A, -J. (92) 

Writing an = a = —022, we now have 



1 



2A 



A = a\^ l^l^-M. (93) 

\|AP-1 ^ ' 



The right-hand side of (|l^) becomes 



2A /,^2 l^ A 



by the Cayley-Hamilton theorem. Thus B must still commute with X. 

We now argue as in Case 1 of Section 5.6 and conclude that the state 
must be one of the slice states (|72|) or (|71D, but with A replaced by 1/A. 

5.8 Non-generic Type 4: TiTa^ = Al. 

The assumption (Gen 4) stated that the only matrices that commute with 
the 2x2 matrix X = TiT2~^ are linear combinations of 1 and X itself. This 
fails only if X is a multiple of the identity, in which case T2 = ATi and the 
state is factorisable: 

I^) = (IT) + A|i))5^ti.,#.)|^,), (95) 



so that particle 1 is a bystander. 
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5.9 Non-generic Type 5 

The assumption (Gen 5) was the statement that al + (3TiT2^^ is not her- 
mitian unless /3 = 0. Suppose this is not true, i.e., 

TiTa"^ = ul + vB (96) 

where u and v are complex scalars and B is hermitian and traceless. To 
analyse states of this form, let us assume that the basis states of particle 1 
have been chosen by means of a Schmidt decomposition of the three-particle 
state l^'), so that the two-particle states 

m = Y.^l^M'^)\i^3) (97) 

and 

\^2) = Y.^2^M^m) (98) 

are orthogonal. Let us also suppose that the basis states of particles 2 and 3 
have been chosen so that T2 is diagonal. Writing 

^-(s:)- '''' 

we then have 

y pzv q\u — rv) J 

and the orthogonality of and |$2) gives 

P^(m + rt>) + g^(u — rt>) = 0. (101) 

Now from ( [25| ) and the following line, the traceless hermitian matrix C 
is given by 

= -T,-'BT,= ( -\_ -P^Vy ^^Q2) 
\^— g pz r J ^ ' 

Since this is hermitian and p and q are real, p^ = q^. Now ( |101| ) gives us 
u = 0, so the state is 

\^)=p\ T)IT)±li)li)) 

+ P^^IT)[r(|T)|T)TU)U))±^| T)li)+^l i)IT)]. 
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We can choose the upper sign (the state with the lower sign is related to it 
by changing the sign of | is)). Then T2 is a multiple of the identity and Ti 
is hermitian, so both T-matrices can be simultaneously diagonalised. Since 
trTi = 0, this gives a state of the form 

1^) = ;^cosa| i) (I T)| T) + I i)| i)) + ^sinad T)| T) - | i)| i)) • 

Relabelling particles 1 and 2 gives 

|vi/) = i=|t)(cosa|i)|T) + Sinai T)|T)) 

+ -^| i)(cosa| i)| i)-sina| t)| i)) 

= ^1 T)l ni T) + ^1 i) (cos2«| ni i) - sin2a| i)) 

where | f) = cosa| J.) + sin a;| t) and | i') = — sin a\ |) + cos a| t)- This is a 
slice ridge state. 

This completes the classification theorem. ■ 

6 A bestiary of atypical pure states of three 
spin-| particles. 

In this section we will summarise the findings of the previous section by de- 
scribing all pure three-particle states with exceptional types of entanglement. 
We will describe their place in the space of all pure three-particle states, us- 
ing the canonical form of Linden, Popescu and Schlienz (henceforth called 



the LPS normal form) from [|T3| These authors pointed out that any 
normalised three-particle state can be brought by local unitary operations to 
the form 

cosal T)(cos/?| T)|T) + sin/?| i)| i)) 
+ sina| i)(-tsin/3| T)| T) + tcos/3| i)| I) + s\ T)| i) + ^| i)| T)) (103) 

where a and (3 are angles lying between and j, t and s are real and 
positive, and 

+ t'^ + \z\^ = 1. (104) 

In accord with our remark at the end of section 3, there are five independent 
parameters (the sixth being the norm which we are taking to be 1). States 
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with different values of these five parameters are locally inequivalent, except 
that when r = or s = we may change the phase of z, which may therefore 
be taken to be real and positive; and when a = all values of (s, t, z) give 
the same state. 

We will also give an indication of the exceptional nature of these states 
and their physical significance by calculating their 2-tangles and 3-tangles. 
These invariants, which were introduced by Wootters 0, quantify how much 
of the entanglement is contained in particular pairs and how much is an 
essential property of the full set of three particles. Formulae for them were 
given by Coffman, Kundu and Wootters [Q. For a pure three-particle state, 
the 2-tangle of particles A and B is 

TAB = [max{Ai - A2 - A3 - A4, 0}]2 (105) 

where (Ai, A2, A3, A4) are, in decreasing order of magnitude, the positive 
square roots of the eigenvalues of 

PabPab = Pab{pab - Pa - + 1), (106) 

Pab being the reduced density matrix of the pair {A,B), obtained from 
|\E')(^| by tracing over particle C, while Pa, Pb are the reduced density ma- 
trices of particles A and B. The 3-tangle is 

tabc = 4 det Pa - tab - tac (107) 

which can be shown [^] to be invariant under permutations of A, B and C . 
The exceptional states are as follows. 

6.1 Bystander States 

These are states which factorise as the product of a one-particle state and 
a two-particle state, so that the one particle is a bystander. They occur 
when the LPS parameters have the values a = 0or/5 = 0, s = t = Oor 
/3 = 0, s = z = 0. The state given by a = 0, namely 

|T)(cos/3| T)|T)+sin/5| DID) 

has the two-dimensional stabiliser 

unless 13 = J when the two-particle state is maximally entangled and the 
stabiliser is four-dimensional: 

(e^'^, f/, V, W) = (e^^ e-^^"^ V,V) 
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or P — 0, when the state is completely factorisable and the stabihser is 
three-dimensional: 

(e*^, U, V, W) = (e*'^, e'^^^ e''^"^ 6^""^) 

with ip + 9 + K, + rj = 0. 

The 2-tangles and 3-tangle of this state are 

Ti2 = Ti3 = 0, T23 = sin^ 2/3, 



Tl23 = 0. 

6.1.1 The General Slice State 

These are states given by 

and their relatives obtainable by permuting the particles: explicitly, 

P\ m) + q\ m)+r\ UT), 

p\ m)+q\ m)+r\ iU), 

p\ m)+q\ m)+r|TU)- 

Such states occur among the LPS normal forms when a ^ and any two of 
/3, s and z are zero. They have one-dimensional stabihsers each consisting of 
four circles; for the first state hsted above, the stabihser contains 

(e'^ U, V, W) = e'^'^^^ £2e-'''^^ £i£2l) • (108) 
where £i, £2 = ±1- Its tangle invariants are 

ri2 = 4|p^r|^ Ti3 = r23 = 0, (109) 

6.2 The Maximal Slice State, or "Slice Ridge" 

These states, which are those slice states that have maximal values of two out 
of the three two-particle von Neumann entropies, occur when a Slice state 
has |pp = |gp + |rp = 1/2, i.e. a = ^ in the Linden-Popescu normal form. 
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In addition to the other shce stabihser elements ( |108| ), they have a further 
one-dimensional set of stabiliser elements given, for states in LPS normal 
form with (3 = 0, s = 0, t = cos 7 and z = sin 7, by 



(e^'^, U, V, W) 



e \ / e \ . /sin 7 cos 7 
e J \—e J ycos7 — sm7 

where ei, = ±1 and 6 can take any value between and 2it. 

The tangles of these states continue to be given by (|109|) and ( |11CI| ). Note 



that for given p, the maximum 3-tangle occurs at r = 0, when the state 
belongs to the following class and the stabiliser becomes two-dimensional. 

6.3 Generalised GHZ States 

Occurring at the boundary of the set of slice states, these states are of the 
form 

p\m) + q\m) (biT^ki). 

They have two-dimensional stabilisers 

(e^'^, U, V, W) = (±1, e^^"^ e^''"^ e'^"^) (111) 

with 6 + K, + ri = 0oTTT. In LPS normal form, these states have f3 = 0, s = 
and z = 0. These states have pure three-particle entanglement, since each of 
their two-particle density matrices is 

P12 = Pi3 = P23 = bri TT)(TT I + IgH I 

which is separable. This is shown by the tangle invariants: 

T12 = Ti3 = T23 = 0, 

ri23=4b|>p. (112) 

6.4 The true GHZ State 

This occupies the same position among the generalised GHZ states as the slice 
ridge states among the general slice states, occurring when \p\ = \q\ [a = j 
in LPS normal form), which maximises the 3-tangle ( p.l2| ). In addition to the 
stabiliser elements (|111|) , it has the further two-dimensional set of stabiliser 



elements 

(e^^, U, V, W) = {±i, icr2e^^"% ia2e^''"^ zaae^''"^) 
with 9 + K, + T] = 0. 
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6.5 The Singular Tetrahedral, or "Beechnut" State 

We call "tetrahedral" states of the form 

^ (o p) ' ^ (r 

since when the eight coefficients tij^ are laid out in a 2 x 2 cubic array, these 
states have zero entries except at the vertices of a tetrahedron. If all four of 
a, b, c, d are non-zero, the state is generic. If one of them is zero, say s = 0, 
the state is of the form 

P\ ni) + q\ m)+r| UT) 
which has the one-dimensional stabiliser 

(e^^, U, V, W) = (e*'^, 6^^'"% 6^^'"% e^^'^^) 
Its tangle invariants are 

2 
2 

7-23 "'^'^ 



bl 


2 






2 




bl 




\r\ 




2 




\q\ 




r ' 



ri23 = 0. 

These states are, in a sense, the opposites of the generalised GHZ states: 
their entanglement is concentrated in two-particle entanglement, and they 
have no three-particle entanglement. 



7 Conclusion 

We have mapped the full range of entanglement properties of pure states of 
three spin-^ particles, using their behaviour under local unitary transforma- 
tions as an indicator. We have identified all the types of exceptional states, 
and have shown that these states will have a special relation to certain local 
invariants. In future work we hope to identify these invariants, and to study 
more fully the variation of known invariants, such as the two-particle von 
Neumann entropies, with respect to entanglement type. 
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8 Appendix: The Bestiary's Family Album 



In this collection of figures we reproduce some graphs of the two-particle 
subsystem von Neumann entropies for the various kinds of non-generic state. 
First of all, let us look at the space of all possible pure states of three spin-| 
particles, a shape we nicknamed "The Pod" in figure 0. Then there are the 
Slice States in figure ^ and the Beechnut states in figure 0. 




Figure 1: The Pod: Here is the space of all possible pure states of 
three spin-j particles, shown from two angles. The "hiccup" or seam in 
the parametrisation lines is not a graphical artefact, it is the line where the 
pod surface ceases to be identical to the beechnut surface (see figure 
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Figure 2: The Slice States: The von Neumann entropies for all three 
sets of slice states. The central spine linking all three fins is the suhspace 
of generalised GHZ states, with the maximally entangled GHZ state at the 
top end, and the spin eigenstate at the bottom. The outside corners are the 
three possible two-particle maximally entangled states ( with the other particle 
a bystander), and the edges running from those corners to the spin eigenstate 
have non-maximal two-particle entanglement, hut are still bystander states. 
The edges that run from the points of each maximal two-particle entanglement 
to the maximal GHZ state are the slice ridges. 




Figure 3: The Beechnut: Here's how the Beechnut states got their name. 
These are the same graph, seen from two angles. Note that the dome at the 
top of the Beechnut doesn't reach the maximally entangled GHZ state. 
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